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Abstract. The Liouville property of a complete Riemannian manifold M 
(i.e., the question whether there exist non-trivial bounded harmonic functions 
on M) attracted a lot of attention. For Cartan-Hadamard manifolds the role 
of lower curvature bounds is still an open problem. We discuss examples of 
Cartan-Hadamard manifolds of unbounded curvature where the limiting angle 
of Brownian motion degenerates to a single point on the sphere at infinity, but 
where nevertheless the space of bounded harmonic functions is as rich as in 
the non-degenerate case. To see the full boundary the point at infinity has to 
be blown up in a non-trivial way. Such examples indicate that the situation 
concerning the famous conjecture of Greene and Wu about existence of non- 
trivial bounded harmonic functions on Cartan-Hadamard manifolds is much 
more complicated than one might have expected. 
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1. Introduction 

The study of harmonic functions on complete Riemannian manifolds, i.e., the 
solutions of the equation Am = where A is the Laplace-Beltrami operator, lies 
at the interface of analysis, geometry and stochastics. Indeed, there is a deep 
interplay between geometry, harmonic function theory, and the long-term behaviour 
of Brownian motion. Negative curvature amplifies the tendency of Brownian motion 
to move away from its starting point and, if topologically possible, to wander out to 
infinity. On the other hand, non-trivial asymptotic properties of Brownian paths for 
large time correspond with non-trivial bounded harmonic functions on the manifold. 

There is plenty of open questions concerning the richness of certain spaces of 
harmonic functions on Riemannian manifolds. Even in the case of simply connected 
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negatively curved Riemannian manifolds basic questions are still open. For instance, 
there is not much known about the following problem posed by Wu in 1983. 

Question 1.1 (cf. Wu [27] p. 139). If M is a simply-connected complete Riemann- 
ian manifold with sectional curvature < — c < 0, do there exist n bounded harmonic 
functions (n — dim M) which give global coordinates on Ml 

It should be remarked that, under the given assumptions, it is even not known 
in general whether there exist non-trivial bounded harmonic functions at all. This 
question is the content of the famous Greene- Wu conjecture which asserts existence 
of non-constant bounded harmonic functions under slightly more precise curvature 
assumptions. 

Conjecture 1.2 (cf. Greene- Wu [10] p. 767). Let M be a simply- connected complete 
Riemannian manifold of non-positive sectional curvature and Xq € M such that 

Sectf < -cr( X y 2 for all x£M\K 

for some K compact, c > and r — dist(xo, •)■ Then M carries non-constant 
bounded harmonic functions. 

Concerning Conjecture 1 1 . 21 substantial progress has been made since the pioneer- 
ing work of Anderson [3], Sullivan [5J, and Anderson-Schoen [4J. Nevertheless, the 
role of lower curvature bounds is far from being understood. 

From a probabilistic point of view, Conjecture 11.21 concerns the eventual be- 
haviour of Brownian motion on these manifolds as time goes to infinity. Indeed, for 
any Riemannian manifold, we have the following probabilistic characterization. 

Lemma 1.3. For a Riemannian manifold (M,g) the following two conditions are 
equivalent: 

i) There exist non-constant bounded harmonic functions on M . 

ii) BM has non-trivial exit sets, i.e., if X is a Brownian motion on M then 
there exist open sets U in the 1-point compactification M of M such that 

F{X t € U eventually} ^ or 1. 

More precisely, Brownian motion X on M may be realized on the space C (R+ , M) 
of continuous paths with values in the 1-point compactification M of M, equipped 
with the standard filtration J?" t = a{X s = pr s \s < t] generated by the coordinate 
projections pr s up to time t. Let £ = sup{£ > : X t G M} be the lifetime of X and 
let J^"inv denote the shift-invariant cr-field on C(R+,M). Then there is a canoni- 
cal isomorphism between the space M\, (M) of bounded harmonic functions on M 
and the set b^ nrv of bounded J?i nv -measurable random variables up to equivalence, 
given as follows: 

(1.1) J&(M) &#W~, « * Hm(u ° X t ). 

(Bounded shift-invariant random variables are considered as equivalent, if they 
agree P x -a.e., for each x € M.) Note that the isomorphism (|1.1|) is well defined by 
the martingale convergence theorem, and that the inverse map to (|1.1[) is given by 
taking expectations: 

(1.2) 6JW~3#i — >ueJf h (M) where u{x) :=E X [H]. 
In particular, 

u(x) := P x {X t S U eventually} 
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is a bounded harmonic function on M, and non-constant if and only if U is a 
non-trivial exit set. 

Now let {M,g) be a Cartan-Hadamard manifold, i.e., a simply connected com- 
plete Riemannian manifold of non-positive sectional curvature. (All manifolds are 
supposed to be connected). In terms of the exponential map exp : T Xo M M 
at a fixed base point xo £ M, we identify p: R™ = T Xo M M. Via pullback 
of the metric on M, we get an isometric isomorphism (M,g) = (W n ,p*g), and in 
particular, M \ {xo} = ]0,oo[ x S n . In terms of such global polar coordinates, 
Brownian motions X on M may be decomposed into their radial and angular part, 

X t = (r(X t ),#(X t )) 

where r(X t ) = dist(a;o,X t ) and where "0(X t ) takes values in S n . 

For a Cartan-Hadamard manifold M of dimension n there is a natural geometric 
boundary, the sphere at infinity Soo(M), such that MU5oo(M) equipped with the 
cone topology is homeomorphic to the unit ball B C R d with boundary dB = S d ^ 1 , 
cf. [9], [5]. In terms of polar coordinates on M, a sequence (r„,^„)„ e N of points in 
M converges to a point of Soo(M) if and only if r n — > oo and i9„ — > $ G S n ~ 1 . 

Given a continuous function / : Soo (M) — > R the Dirichlet problem at infinity is 
to find a harmonic function /i : M — > R which extends continuously to (M) and 
there coincides with the given function /, i.e., 

h\S oa {M) = f. 

The Dirichlet problem at infinity is called solvable if this is possible for every such 
function /. In this case a rich class of non-trivial bounded harmonic functions on 
M can be constructed via solutions of the Dirichlet problem at infinity. 

In 1983, Anderson |3] proved that the Dirichlet problem at infinity is indeed 
uniquely solvable for Cartan-Hadamard manifolds of pinched negative curvature, 
i.e. for complete simply connected Riemannian manifolds M whose sectional cur- 
vatures satisfy 

-a 2 < Sectf < -b 2 for all x E M, 

where a 2 > b 2 > arc arbitrary constants. The proof uses barrier functions and 
Perron's classical method to construct harmonic functions. Along the same ideas 
Choi [7] showed in 1984 that in rotational symmetric case of a model (M, g) the 
Dirichlet problem at infinity is solvable if the radial curvature is bounded from 
above by — A/(r 2 log(r)). Hereby a Riemannian manifold (M,g) is called model if 
it possesses a pole p S M and every linear isometry (p: T p M — > T p M can be realized 
as the differential of an isometry M — > M with Q(p) = p. Choi [7] furthermore 
provides a criterion, the convex conic neighbourhood condition, which is sufficient 
for solvability of the Dirichlet problem at infinity. 

Definition 1.4 (cf. Choi [7]). A Cartan-Hadamard manifold M satisfies the convex 
conic neighbourhood condition at x £ 6*00 (M) if for any y 6 Soo(Al), y ^ x, there 
exist subsets V x and V y C M U S 00 (M) containing x and y respectively, such that 
V x and V y are disjoint open sets of M U Soo(M) in terms of the cone topology and 
V X C\M is convex with C 2 -boundary. If this condition is satisfied for all x £ 5oo(M), 
M is said to satisfy the convex conic neighbourhood condition. 

It is shown in Choi [7] that the Dirichlet problem at infinity is solvable for a 
Cartan-Hadamard manifold M with sectional curvature bounded from above by 
— c 2 for some c > 0, if M satisfies the convex conic neighbourhood condition. 
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In probabilistic terms, if the Dirichlet problem at infinity for M is solvable and 
almost surely 

X c := \imX t 

exists in Soo(M), where (Xt)t<( is a Brownian motion on M with lifetime £, the 
unique solution h: M — > R to the Dirichlet problem at infinity with boundary 
function / is given as 

(1.3) h(x)=E[foX ((x) (x)]. 

Here X(x) is a Brownian motion starting at x G M. 

Conversely, supposing that for Brownian motion X{x) on M almost surely 

lim XAx) 

exists in Sqc^M^) for each x tE one may consider the harTnojiic jneasure fi x on 
5oo(M), where for a Borel set U C S^M) 

(1.4) (i x (U) :=F{X ({x) (x) e U). 
Then, for any Borel set U C Soo(M), the assignment 

x i ^ Hx(U) 

defines a bounded harmonic function hu on M. By the maximum principle hjj 
is either identically equal to or 1 or takes its values in ]0, 1[. Furthermore, all 
harmonic measures [i x on Sqc^M) are equivalent. Thus, by showing that the har- 
monic measure class on Soo (M) is non-trivial, we can construct non-trivial bounded 
harmonic functions on M. To show that, for a given continuous boundary function 
/: Soo(M) — > K, the harmonic function 

(1.5) h(x) = f f(y)ti x (dy). 

extends continuously to the boundary Soo(M) and takes there the prescribed bound- 
ary values /, we have to show that, whenever a sequence of points Xi in M converges 
to Xoo € Soo(M) then the measures fx Xi converge weakly to the Dirac measure at x^. 
In particular, for a continuous boundary function f\S OQ (M), the unique solution to 
the Dirichlet problem at infinity is given by formula f) 1 . 5ft . 

The first results in this direction have been obtained by Prat [231 US] between 
1971 and 1975. He proved that on a Cartan-Hadamard manifold with sectional 
curvature bounded from above by a negative constant — fc 2 , k > 0, Brownian motion 
is transient, i.e., almost surely all paths of the Brownian motion exit from M at 
the sphere at infinity [25] . If in addition the sectional curvatures are bounded from 
below by a constant — K 2 , K > k, he showed that the angular part $(X t ) of the 
Brownian motion almost surely converges as t — > Q. 

In 1976, Kifer [17] presented a stochastic proof that on Cartan-Hadamard mani- 
folds with sectional curvature pinched between two strictly negative constants and 
satisfying a certain additional technical condition, the Dirichlet problem at infinity 
can be uniquely solved. The proof there was given in explicit terms for the two 
dimensional case. The case of a Cartan-Hadamard manifold (M, g) with pinched 
negative curvature without additional conditions and arbitrary dimension was fi- 
nally treated in Kifer [18] . 

Independently of Anderson, in 1983, Sullivan [26] gave a stochastic proof of the 
fact that on a Cartan-Hadamard manifold with pinched negative curvature the 
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Dirichlet problem at infinity is uniquely solvable. The crucial point has been to 
prove that the harmonic measure class is non-trivial in this case. 

Theorem 1.5 (Sullivan [26j). The harmonic measure class on Soo(M) — d(M U 
Soo(M)) is positive on each non-void open set. In fact, if Xi in M converges to Xoo 
in Soo(M), then the Poisson hitting measures fj, Xi tend weakly to the Dirac mass 
at • 

In the special case of a Riemannian surface M of negative curvature bounded 
from above by a negative constant, Kendall [Hj gave a simple stochastic proof that 
the Dirichlet problem at infinity is uniquely solvable. He thereby used the fact 
that every geodesic on the Riemannian surface "joining" two different points on 
the sphere at infinity divides the surface into two disjoint half-parts. Starting from 
a point x on M, with non-trivial probability Brownian motion will eventually stay 
in one of the two half-parts up to its lifetime. As this is valid for every geodesic and 
every starting point x, the non-triviality of the harmonic measure class on S^M) 
follows. 

Concerning the case of Cartan-Hadamard manifolds of arbitrary dimension sev- 
eral results have been published how the pinched curvature assumption can be 
relaxed such that still the Dirichlet problem at infinity for M is solvable, e.g. [T3] 
and [12]. 

Theorem 1.6 (Hsu |12j). Let (M,g) be a Cartan-Hadamard manifold. The Dirich- 
let problem at infinity for M is solvable if one of the following conditions is satisfied: 

(1) There exists a positive constant a and a positive and nonincreasing function 
h with J Q rh(r) dr < oo such that 

-h{r{x)) 2 e 2ar{x) < Ricf and Sectf < -a 2 for all x e M. 

(2) There exist positive constants r$, a > 2 and j3 < a — 2 such that 

~r(x) 2fj < Ricf and Sectf < - "^"^ 

x ~ r(x) 2 

for all x £ M with r(x) > rg. 

It was an open problem for some time whether the existence of a strictly negative 
upper bound for the sectional curvature could already be a sufficient condition for 
the solvability of the Dirichlet problem at infinity as it is true in dimension 2. 
In 1994 however, Ancona [5] constructed a Riemannian manifold with sectional 
curvatures bounded above by a negative constant such that the Dirichlet problem 
at infinity for M is not solvable. For this manifold Ancona discussed the asymptotic 
behaviour of Brownian motion. In particular, he showed that Brownian motion 
almost surely exits from M at a single point com on the sphere at infinity. Ancona 
did not deal with the question whether M carries non-trivial bounded harmonic 
functions. 

Borbely [5] gave another example of a Cartan-Hadamard manifold with curvature 
bounded above by a strictly negative constant, on which the Dirichlet problem at 
infinity is not solvable. Borbely does not discuss Brownian motion on this manifold, 
but he shows using analytic methods that his manifold supports non-trivial bounded 
harmonic functions. 

This paper aims to give a detailed analysis of manifolds of this type and to 
answer several questions. It turns out that the manifolds of Ancona and Borbely 
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share quite similar properties, at least from the probabilistic point of view. On 
both manifolds the angular behaviour of Brownian motion degenerates to a single 
point, as Brownian motion drifts to infinity. In particular, the Dirichlet problem 
at infinity is not solvable. Nevertheless both manifolds possess a wealth of non- 
trivial bounded harmonic functions, which come however from completely different 
reasons than in the pinched curvature case. Since Borbely's manifold is technical 
easier to handle, we restrict our discussion to this case. It should however be noted 
that all essential features can also be found in Ancona's example. 

Unlike Borbely who used methods of partial differential equations to prove that 
his manifold provides an example of a non-Liouville manifold for that the Dirichlet 
problem at infinity is not solvable, we are interested in a complete stochastic de- 
scription of the considered manifold. In this paper we give a full description of the 
Poisson boundary of Borbely's manifold by characterizing all shift-invariant events 
for Brownian motion. The manifold is of dimension 3, and the er-algebra of invari- 
ant events is generated by two random variables. It turns out that, in order to see 
the full Poisson boundary, the attracting point at infinity to which all Brownian 
paths converge needs to be "unfolded" into the 2-dimensional space Ex S 1 . 

The asymptotic behaviour of Brownian motion on this manifold is in sharp con- 
trast to the case of a Cartan-Hadamard manifold of pinched negative curvature. 
Recall that in the pinched curvature case the angular part $(X) of X carries all rel- 
evant information and the limiting angle lim^oo $(X t ) generates the shift-invariant 
a- field of X. All non-trivial information to distinguish Brownian paths for large 
times is given in terms of the angular projection of X onto Soo(M), where only the 
limiting angle matters. As a consequence, any bounded harmonic function comes 
from a solution of the Dirichlet problem at infinity, which means on the other hand 
that any bounded harmonic function on M has a continuous continuation to the 
boundary at infinity. 

More precisely, denoting by h(M) the Banach space of bounded harmonic func- 
tions on a Cartan-Hadamard manifold M, we have the following result of Anderson. 

Theorem 1.7 ([5]). Let (M,g) be a Cartan-Hadamard manifold of dimension d, 
whose sectional curvatures satisfy —a 2 < Sect^ < — b 2 for all x e Af. Then the 
linear mapping 

P: L oo (S 0O (M),n)^h(M), 
(L6) f^P(f), P(f)(x):= [ fdfi x 

is a norm-nonincreasing isomorphism onto h(M). 

In the situation of Borbely's manifold Brownian motion almost surely exits from 
the manifold at a single point of the sphere at infinity, independent of its starting 
point x at time 0. In particular, all harmonic measures fi x are trivial, and harmonic 
functions of the form P(f) are necessarily constant. On the other hand, we are 
going to show that this manifold supports a huge variety of non-trivial bounded 
harmonic functions, necessarily without continuation to the sphere at infinity. It 
turns out that the richness of harmonic functions is related to two different kinds 
of non-trivial exit sets for the Brownian motion: 

1. Non-trivial shift-invariant information is given in terms of the direction along 
which Brownian motion approaches the point at the sphere at infinity. Even if 
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there is no contribution from the limiting angle itself, angular sectors about the 
limit point allow to distinguish Brownian paths for large times. 

2. As Brownian motion converges to a single point, in appropriate coordinates 
the fluctuative components (martingale parts) of Brownian motion for large times 
are relatively small compared to its drift components. Roughly speaking, this means 
that for large times Brownian motion follows the integral curves of the deterministic 
dynamical system which is given by neglecting the martingale components in the 
describing stochastic differential equations. We show that this idea can be made 
precise by constructing a deterministic vector field on the manifold M such that 
the corresponding flow induces a foliation of M and such that Brownian motion 
exits M asymptotically along the leaves of this foliation. 

To determine the Poisson boundary the main problem is then to show that 
there are no other non-trivial exit sets. This turns out to be the most difficult 
point and will be done by purely probabilistic arguments using time reversal of 
Brownian motion. To this end, the time-reversed Brownian motion starting on the 
exit boundary and running backwards in time is investigated. 

It is interesting to note that on the constructed manifold, despite of diverging 
curvature, the harmonic measure has a density (Poisson kernel) with respect to the 
Lebesgue measure. Recall that in the pinched curvature case the harmonic measure 
may well be singular with respect to the surface measure on the sphere at infinity 
(see [17], [22], [15], [8] for results in this direction). Typically it is the fluctuation 
of the geometry at infinity which prevents harmonic measure from being absolutely 
continuous. Pinched curvature alone does in general not allow to control the angular 
derivative of the Riemannian metric, when written in polar coordinates. 

The paper is organized as follows. In Section [2] we give the construction of our 
manifold M which up to minor technical modifications is the same as in Borbely [6] : 
we define M as the warped product 

M := (HDL) XgS 1 , 

where L is a unit-speed geodesic in the hyperbolic space H 2 of constant sectional 
curvature —1 and H is one component of H 2 \ L. The Riemannian metric 7 on 
M is the warped product metric of the hyperbolic metric on H coupled with the 
(induced) Euclidean metric on S 1 via the function g: H U L — ► K + , 

ds 2 M = ds^ 2 + g ■ ds 2 sl . 

By identifying points (i, a\) and (£, 0:2) with I S L and a\, a.2 € S 1 and choosing the 
metric "near" L equal to the hyperbolic metric of the three dimensional hyperbolic 
space H 3 the manifold M becomes complete, simply connected and rotationally 
symmetric with respect to the axis L. We specify conditions the function g has 
to satisfy in order to provide a complete Riemannian metric on M for which the 
sectional curvatures are bounded from above by a negative constant, and such that 
the Dirichlet problem at infinity is not solvable. As the construction of the function 
g is described in detail in [5] we mainly sketch the ideas and refer to Borbely 
for detailed proofs. We explain which properties of g influence the asymptotic 
behaviour of the Brownian paths. 

The probabilistic consideration of the manifold M starts in Section [3] We spec- 
ify the defining stochastic differential equations for the Brownian motion X on 
M, where we use the component processes R, S and A of X with respect to the 
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global coordinate system {(r,s,a)\r € R+, s <E R, a € [0, 2tt)} for M. The non- 
solvability of the Dirichlet problem at infinity is then an immediate consequence of 
the asymptotic behaviour of the Brownian motion (see Theorem I3.2[) . 

It is obvious that the asymptotic behaviour of the Brownian motion on M is the 
same as in the case of the manifold of Ancona. In particular, it turns out that the 
component R of the Brownian motion X almost surely tends to infinity as t — > £. It 
is a remarkable fact that in contrast to the manifold of Ancona, where the Brownian 
motion almost surely has infinite lifetime, we can show that on M the lifetime C of 
the Brownian motion is almost surely finite. 

In Section [4] we start with the construction of non-trivial shift-invariant events. 
To this end, we consider a time change of the Brownian motion such that the drift of 
the component process R of X equals t, i.e. such that the time changed component 
R behaves then comparable to the deterministic curve R+ — > R+, t i— > ro + 1. We 
further show that for a certain function q, the process 

R t 



Z t := S t - / q(r) dr 
Jo 



almost surely converges in R when t — > £. 

It turns out that the non-trivial shift-invariant random variable Aq := lim t _>£ A t 
can be interpreted as one-dimensional angle which indicates from which direc- 
tion the projection of the Brownian path onto the sphere at infinity attains the 
point L(po). The shift-invariant random variable Z^ indicates along which surface 
of rotation C SQ x S 1 inside of M the Brownian paths eventually exit the manifold M. 
Thereby C So is the trajectory starting in (0, s ) € R+ x R of the vector field 

V:=Uq{r)° 
or os 

Section [5] finally gives a description of the Poisson boundary. The main work 
in this part is to verify that there are no other invariant events than the specified 
ones. This is achieved by means of arguments relying on the time-inversed process. 

As a consequence a complete description of all non-trivial bounded harmonic 
functions on the manifold M is obtained. Concerning "boundary Harnack inequal- 
ity" the class of non-constant bounded harmonic functions shares an amazing fea- 
ture: on any neighbourhood in M of the distinguished point at the horizon at 
infinity, a non-trivial bounded harmonic function on M attains each value lying 
strictly between the global minimum and the global maximum of the function. 

Finally it is worth noting that our manifold also provides an example where the 
(T-field of terminal events is strictly bigger than the cr-field of shift-invariant events. 



2. Construction of a CH manifold with a sink of curvature at 

INFINITY 

Let L C H 2 be a fixed unit speed geodesic in the hyperbolic half plane 

H 2 = {( a; ,y)eR 2 |y>0} 

equipped with the hyperbolic metric ds^p of constant curvature —1. For our pur- 
poses one can assume without loss of generality L := {(0, y)\y > 0} to be the 
positive y-axis. Let H denote one component of H 2 \ L and define a Riemannian 
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manifold M as the warped product: 

M := (H U L) XgS 1 , 

with Riemannian metric 

ds 2 M = cLsjjp + g ■ ds 2 S i , 

where g : H U L — > R + is a positive C°° -function to be determined later. By 
identifying points (£, oti) and (£,012) with I € L and ax,Qf2 ^ ""1 we ma ke M a 
simply connected space. 

We introduce a system of coordinates (r, s, a) on M, where for a point p € M 
the coordinate 7' is the hyperbolic distance between p and the geodesic L, i.e. the 
hyperbolic length of the perpendicular on L through p. The coordinate s is the 
parameter on the geodesic {L(s) : s 6 ]— 00, oo[}, i.e. the length of the geodesic 
segment on L joining L(0) and the orthogonal projection L(s) of p onto L. Further- 
more, a £ [0, 2ir[ is the parameter on S 1 when using the parameterization e la . We 
sometimes take a £ R, in particular when considering components of the Brownian 
motion, thinking of R as the universal covering of S 1 . 

In the coordinates (r, s, a) the Riemannian metric on M \ L takes the form 

(2.1) ds 2 M = dr 2 + h(r)ds 2 + g(r, s)da 2 

where h(r) = cosh 2 (r). 

Let g(r,s) := sinh (r) for r < 1/10 (the complete definition is given below), 
then the above metric smoothly extends to the whole manifold M, where M is now 
rotationally symmetric with respect to the axis L and for r < 1/10 isometric to 
the three dimensional hyperbolic space H 3 with constant sectional curvature — 1, 
cf. [6]. From that it is clear that the Riemannian manifold (M,g) is complete. 

2.1. Computation of the sectional curvature. From now on we fix the basis 
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for the tangent space T P M in p £ M. Herein the Christoffel symbols of the Levi- 
Civita connection can be computed as follows - the indices refer to the correspond- 
ing tangent vectors of the basis: 

1 h' n' 

pl _ x j / y~2 _ r 2 _ n r r 3 _ r 3 _ i/r 

122 rl 12 21 ~ 2/? il3_i31 ~2^' 



r 1 — n' r 2 — —Q' s r 3 _ r 3 _ d' s 

i 3 3 - 2 9r, 133- 2h > 1 2 3 - i 32- 2ff > 

all others equal 0. Herein g' r denotes the partial derivative of the function g with 
respect to the variable r, g' s the partial derivative with respect to s, etc. 
For the computation of the sectional curvatures Sect M write X = (a;i, #2, £3) and 
Y = (2/1)2/2,2/3) for tangent vectors X, Y e T p M in terms of the basis 81,82,83. 
Then one gets 

(R(X, Y)Y,X) = (- l -K r + {xiV2 x 2Vl f 



-.A 



= :B 

±g>J a -lg>X + \^)(.r. 2 ,,,- r,,, 2 f 



:C 



+ ■2- [ ^ I -^IAl | iMi ) - j (.;■,;,,. -.r,// 2 ). 



:D 



as well as 

\\X A Y\\ 2 = h(xiy 2 - x 2 y x f + g{xiy 3 - x 3 yi) 2 + gh(x 2 y 3 - x 3 y 2 ) 2 . 
We conclude that the manifold M has strictly negative sectional curvature, i.e. 



-k 2 > Sect M (span{X,r}) = 



X AY 



for some k > 0, all X, Y € T p M and all p <G M, if and only if the following 
inequalities hold: 

(2.2) h'; r - ]&. > k x 



2' rr Ah 
(2-3) ^_^>fc» fl> 

(2-4) \g'L + \g'X~Y^>k 2 gh, 

2 4 4 5 



1 ( I lg_X lg^\ 2 (g'; r (g' r ) 2 

Mrs ' A U A I — 



g 2 h V ^ Ah A g J ~ \ 2g Ag 2 

(2-5) 



(g'L g'rK {g's? , 

^2 5 /i Agh Ag 2 h 
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FIGURE 2. Sphere at Infinity S^M) 



This is explained by the fact that the quadratic form 

q(X, Y, Z) := (A + k 2 h)X 2 + (B + k 2 g)Y 2 + (C + k 2 gh)Z 2 + 2DYZ 
is non-positive for all X, Y, Z £ R if and only if 
-A > k 2 h and - B > k 2 g and - C > k 2 gh and D 2 < (B + k 2 g)(C + k 2 gh). 

2.2. The sphere at infinity S oa (M). As it is obvious that, for (s, a) € Rx [0, 2tt[ 
fixed, the curves 

7 SQ : M + -> M, r i-> (r, s, a) 

form a foliation of M by geodesic rays, we can describe the sphere at infinity Soo(M) 
as the union of the "endpoints" (i.e. equivalence classes) r y sa (+oo) of all geodesic 
rays 7 SQ foliating M together with the equivalence classes L(+oo) and L{— oo) 
determined by the geodesic axis L of M, which sums up to: 

Soo(M) = L(+oo) U { 7sct (+oo) | (s,a) £ K x [ 0, 2vr[} U L(-oo). 

This also explains why it suffices to show that the s-component St of the Brow- 
nian motion X t converges to +oo as t approaches the lifetime £ of the Brownian 
motion, if we want to show that X t converges for t — > Q to the single point L(+oo) 
on the sphere S aD (M) at infinity. 

2.3. Properties of the function g. We give a brief description of the properties 
which the warped product function g : R + xl^ R + in (12. ip has to satisfy to provide 
an example of a Riemannian manifold where the Dirichlct problem at infinity is not 
solvable whereas there exist non-trivial bounded harmonic functions. We give the 
idea how to construct such a function g and refer to Borbely [6 for further details. 

Lemma 2.1. There is a C°° -function g: M+ x R — > R + . (r, s) g(r, s) satisfying 
the following properties: 

(1) g' r > and g' s > 0. 

(2) g(r, s) = sinh 2 (r) for r < 1/10, and for r > 1/10 it holds that: 

,~> h , A 1 // 1 (9r) 2 . K9'r 

g r > h r g and -g rr - > 
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(3) Denoting p(r, s) := g' s /{g' r h) one has {ph)' s > and for all s € 

p(r, s) dr = oo. 



o 



(4) The function p(r, s) satisfies: 

1 / 1 i/i ^ . / . 9 h r 

P - 7777777' Ps ^ 7777777' br I < 7777777 and PPr^ < 



1000 s ~ 1000 - 1000 r 1000 

(5) We further need that 

{ P h)' r > o, ( P hy; r > o 

and that there exists e € ]0, l/4[ smc/i t/ia£ 

(2.6) (ph)'M > p{r Q £ )fe(r) /or r > 2. 

(6) The function p{r, s) finally has to fulfill: 

~i — TTTT = oo 

r p{r,s)h(r) 

for every r > and every s€l. 
The following remark summerizes briefly the significance of these conditions. 

Remark 2.2 (Comments on the properties listed above). 

1. As mentioned in the definition of the Riemannian metric on M , property ([2]) 
assures that the Riemannian metric smoothly extends to the geodesic L and 
that curvature condition (|2.3|) is satisfied for a suitable k. 

2. The conditions in (j4|) stated for p(r, s) assure validity of the curvature 
condition ((23)) - at least for k = 1/1000. 

3. The integral condition ([3]) forces Brownian motion X — (R,S,A) on M 
to converge to the single point L(oo) £ Soo(M) as t — > £ which as an 
immediate consequence implies non-solvability of the Dirichlet problem at 
infinity for M. The given condition assures that the drift term in the 
stochastic differential equation for the component S of the Brownian motion 
X compared with the drift term appearing in the defining equation for the 
component R grows "fast enough" to ensure St going to oo as t approaches 
the lifetime of the Brownian motion. 

In [6], Lemma 2.1, Borbely uses this condition to show that the convex 
hull of any neighbourhood of L(oo) is the whole manifold M. This is a 
natural first step in the construction of a manifold for which the Dirichlet 
problem at infinity is not solvable (cf. the Introduction and [7j, [6]). 

4. Property ([6]) is needed for the Brownian motion to almost surely enter a 
region where St has nonzero and positive drift (see the construction of I 
and S> below) ; otherwise St might converge in K with a positive probability. 

2.4. Construction of the function g. The idea to construct a function g with the 
wanted properties is as follows. As g is given as solution of the partial differential 
equation 

(2.7) g' s (r,s)=p(r,s)h{r)g' r (r,s) 

one has to find an appropriate function p(r, s) and the required initial conditions 
for g to obtain the desired function. 
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We will see later that the construction of the metric on M is similar to that given 
in Ancona [2], as Borbely also defines the function p "stripe- wise" to control the 
requirements for p, g respectively, on each region of the form [r^, rj] x R. Yet he is 
mainly concerned with the definition of the "drift ratio" p, what makes it harder to 
track the behaviour of the metric function g and to possibly modify his construction 
for other situations, whereas Ancona gives a more or less direct way to construct the 
coupling function in the warped product. As a consequence this allows to extend 
his example to higher dimensions and to adapt it to other situations as well. 

We start with a brief description how to construct the function 

Po(r) :=p(r,0), 

as given in [6]: po is defined inductively on intervals [r n , r„+i], where r n+ i — r n > 3 
and n > 3 sufficiently large, see below. We let po :— on [0, 2] and define it on [2, 3] 
as a slowly increasing function satisfying conditions (HJ and (JSJ). For re [3,ri] let 
Po(r) := po(3) be constant, where r\ is chosen big enough such that (j>oh)(ri) > 1 
and r\/h{r\) < 1/1000. On the interval [ri,oo[ we choose the function p to 
be decreasing with limr^oo po(r) — 0, whereas poh is still increasing and strictly 
convex; see [6], Lemma 2.3 and Lemma 2.4. 

On the interval [ri,^] for ri big enough as given below (and in general on 
intervals of the form [f2n— i)^2n]) one extends poh via a solution of the differential 
equation 

y" = V(2y). 

Note that carefully smoothing the function po on the interval [r\ , r± + 1] (and on 
[f2n-i, ?'2n-i + 1] respectively) to become C°° does not interfere with the properties 
(|4]) and iJSJ of po and can be done such that (poh)" > l/(Apoh) is still valid. 

Lemma 2.4 in [6] shows that in fact po = p$h/h decreases on [7"2n-ij ^n] an d 
again by [6], Lemma 2.3, for given ran— l ° ne can choose the upper interval bound 
r2n such that 

r 2 " i 

/ — - dr > 1 for all n, 

Jr 2n -i Poh 

which guarantees property (JH) for po . 

On the interval [r2, r^] with r% sufficiently large as given below and in general on 
intervals of the form [r2„, r , 2„+i], let po — c„/r for some well chosen constant c n . 
This differs from the construction of [5] p. 228 where po is chosen to be constant 
in these intervals, but it does not change the properties of the manifold. As above, 
smoothing on intervals [?*2 n — l,r2 n ] preserves the conditions ([4]), ([5]) and (poh)" > 
s/(poh) for < e < 1/4 small enough and independent of n, but depending on the 
choice of po on [2, 3[. 

If one chooses for given r2, r2„ respectively, the interval bounds r^, T2 n +i re- 
spectively, large enough we can achieve that 

T2n+1 

Po(r) dr > 1 for all n, 

2n 

which finally assures property (jHJ) for p . 

Following Borbely [6], one defines p(r, s) via po(r) as 

p(r,s) := x(r,s)p (r) 

by using a "cut-off function" x( r , s )- Herein x(r, s) is given as x( t j s ) '■= C( s + ^( r ))) 
where ^ is smooth and increasing with £(y) = for y < 0, £(y) = 1/2 for y > 4 and 
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f(fo,8 ) = (l(ra,s ,-£{f(r ,s ))) 



7(r ,s ,0) 



domain & 



graph of s i-» y(r , s ,s) 
integral curve of — Z = —(ph)J^ 



a_ 

Ox 



—£{ro) s o s axis 

Figure 3. Graph of s h-> 7(7*0, s , s) 

< 1/2, £," + (,> 0. The function I satisfies £(0) = 0, l'(r) = on the interval 
[0,2] and t'{r) = e/(poh)(r) on the interval [3, oo[, with the same e as in (12. 6|) . 
Then the two pieces are connected smoothly such that for r > 

.(Pofc)'r(r) 



(2.8) 



f (r) > -E- 



and < £' < . 

Pah 



poh(r 

The function £ is nondecreasing such that lirtv^oo £(r) = oo (see ([5]) in Lemma |2. II) . 
and one can check that p(r, s) satisfies the required properties listed in Lemma l2.ll 
(see [5], p. 228ff). 

For an appropriate initial condition to solve the partial differential equation 
(|2.7p . we set go(r) :— sinh 2 (r) on [0, j^]. On [i,oo[ let <7o(r) be the solution of the 
differential equation 

i . , / n . . > / j_ 

10. 



/' 



- h'f with initial condition f . 
sinh 2 (l/10) r V 10 



( — \ = sinh 2 



Smoothing go yields a C°°-function go such that go = go on [0, ^ — S] U [ 
for an appropriate 5. In particular, there exist d\, d 2 > such that 

1 



(5, oo [ 



(2.9) 



go (r) = d 2 e dl 



d-i sinh 



for all r > 



10 



5. 



The function go serves as initial condition to solve the partial differential equation 

9s( r i s ) =P( r ,s)h(r)g' r (r,s). 



More precisely, it is proven in [6] that all trajectories of the vector field — Z 
ph)J^ + intersect the domain 

9 := {(s, r) G R x E + , s < -£(r)} . 



Letting s i— » (7(7-0, so, s)> s) be the integral curve of — Z satisfying 7(70, so, so) = ro, 
we define /(ro, so) by 



(2.10) 



/(ro, s ) = 7(r , so, -£(f(ro, s ))). 
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Then (see Sect. [5] for details) 

(2-11) g{r , s ) = 9o(f(r , s )). 

3. Asymptotic behaviour of Brownian motion on M 

Let (f2; J^";P) be a filtered probability space satisfying the usual conditions and 
X a Brownian motion on M considered as a diffusion process with generator \^-m 
taking values in the Alexandroff compactification M := M U {oo} of M. Further 
let £ denote the lifetime of X, i.e. Xt{u>) = oo for t > £(w), if £( w ) < °°- 

As we have fixed the coordinate system 

M = {(r,s,a) : r E R+, set, a E [0,2tt[} 

for our manifold M, we consider the Brownian motion X in the chosen coordinates 

as well and denote by {Rt)t<C, {St)t<( and (At)t«; the component processes of 

JL JL JL 

dr 5 ds ' da 



(Xt)t«z- The generator |Am of A is then written in terms of the basis — — 



of TM as: 

K ' } 2 M 2 v 9r 2 /i 9s 2 .g da 2 J \Ah Ag J dr Agh ds' 

Therefore we can write down a system of stochastic differential equations for the 
components R, S and A of our given Brownian motion: 

^qica j D ( Ki R t) g' r (R t ,S t )\ x 

(3 - 13) dRt= {^) + W^)r + ^ 

(3.14) dS t = - ; (ft + -=^= dt^ 2 

ig(Rt,S t )h(R t ) y/hjRt) 

(3.15) dA 4 = 1 

y/g(Rt,S t ) 

with a three-dimensional Euclidean Brownian motion = (W 1 , VF 2 , V^ 3 ). 

As already mentioned, we are going to read the component ^4 of the Brownian 
motion with values in the universal covering R of S 1 . 

Remark 3.1. Inspecting the defining stochastic differential equations for the com- 
ponents R, S and A of the Brownian motion X on M, it is interesting to note 
that the behaviour of the component A t does not influence the behaviour of the 
components Rt and St- The first two equations can be solved independently of the 
third one; their solution then defines the third component of the Brownian motion. 
Hence it is clear that the lifetime C of X t does not depend on the component A t - 
in particular does not depend on the starting point A of A t . We are going to use 
this fact later to prove existence of non-trivial bounded harmonic functions on M. 

As we we are going to see in Section |4*U1 the "drift ratio" p(r,s) = g' s /(g' r h) 
influences the interplay of the components St and Rt of the Brownian motion and 
therefore determines the behaviour of the Brownian paths. For this reason it is 
more convenient to work with a time-changed version Xt of our Brownian motion, 
where the drift of the component R t is just t and the drift of St is essentially given 
by p. This can be realized with a time change (r t ) defined as follows: 

Let 
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and n := T' 1 ^) = inf{s e M+ : T(s) > t} for t < T((). The components R t , S t ,A t 
of the time-changed Brownian motion X t := X Tt are then given for t < £ :— 
by the following system of stochastic differential equations: 

(3.16) dR t = dt + 1 dW 1 

(3.17) dS t = 9's(Rt,S t ) = ^ + 1 dW 2 

(gh' r + g' r h)(R t ,S t ) ^_ {Rt) + 

(3.18) cL4 t = 1 dVF 3 . 

/ 5 (i? t ,S t ) (§(j? t ) + g(i? t ,S t )) 



We are now going to state the main theorem of this chapter which shows that 
from the stochastic point of view the Riemannian manifold (M, 7) constructed by 
Borbely [5] has essentially the same properties as the manifold of Ancona in [2]. 
We further give a stochastic construction of non-trivial bounded harmonic functions 
on M, which is more transparent than the existence proof of Borbely relying on 
Perron's principle. 

Theorem 3.2 (Behaviour of Brownian motion on M). 

(1) For the Brownian motion X on the Riemannian manifold (M, 7) con- 
structed above the following statement almost surely holds: 

lim X t = L(+oo), 

independently of the starting point Xq. In particular the Dirichlet problem 
at infinity for M is not solvable. 

(2) The component A t of Brownian motion X t almost surely converges to a 
random variable which possesses a positive density on S 1 . 

(3) The lifetime ( of the process X is a.s. finite. 

Proof. From Eqs. (|3 . 16[> . (|3.17[) and (|3.18[) we immediately see that the derivatives 
of the brackets of R and S are bounded, and that the drift of S takes its values 
in [0,maxp] and is therefore bounded. As a consequence the process (R, S,A) has 
infinite lifetime. 

Moreover for all [3 > 0, we have almost surely \R t — t\ < (3t eventually. From 
this and the fact that g' r /g > h' r for r > 1/10 we deduce that the martingale parts 
of R and S together with the process A converge as t tends to infinity. 

Next we see from 

9s = P k 9r < 

gh' r + g' r h gh' r + g' r h 

that for t sufficiently large the drift of S is larger that po((l + f3)t) /2. Consequently 
S tends to infinity as t — > 00. To prove ([T]) it is sufficient to establish 

lim X t = L(+oo) 

t — >OG 

and this is a consequence of the fact that St converges to infinity. The proof 
of (J2J is a direct consequence of the convergence of A t to a random variable which 
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conditioned to R and S is Gaussian (in K) and non-degenerate. We are left to 
prove ([3]). Changing back the time of the process X, we get 

(3.19) C= rivt + Tz) dt - 



Ah Ag 
But, for large r and s > 0, we have 

— + — I (r, s) > — > — cosh r sinh r. 
Ah Ag J Ag ~ 2 

As a.s. Rt — t > —fit eventually and as St converges to infinity as t tends to infinity, 
we get the result. □ 

4. NON-TRIVIAL SHIFT-INVARIANT EVENTS 

As explained in the Introduction there is a one-to-one correspondence between 
the er-field J^inv of shift-invariant events for X up to equivalence and the set of 
bounded harmonic functions on M. We are going to use this fact and give a 
probabilistic proof for the existence of non-trivial shift-invariant random variables, 
which in turn yields non-trivial bounded harmonic functions on M. In addition, we 
get a stochastic representation of the constructed harmonic functions as "solutions 
of a modified Dirichlet problem at infinity". In contrast to the usual Dirichlet 
problem at infinity however, the boundary function does not live on the geometric 
horizon S OD (M) = S d ~ 1 (M), and the harmonic functions are not representable in 
terms of the limiting angle of Brownian motion. 

4.1. Shift-invariant variables. In the discussed example it turns out that the 
shift-invariant random variable Aq is non-trivial and can be interpreted as "1- 
dimensional angle" on the sphere S' 00 (M) at infinity. The variable gives the 
direction on the horizon, wherefrom the Brownian motion X converges to the lim- 
iting point L(oo). Despite the fact that the limit X^ 6 Soo(M) itself is trivial, 
and hence the limiting angle € S d ~ 1 (M) as well, Brownian paths for large 

times can be distinguished by their projection onto Sao(M). Taking into account 
that in the pinched curvature case the angular part carries all shift-variant infor- 
mation, one might conjecture that the random variable A^ already generates the 
shift-invariant cr-field J^inv In turn this would imply a stochastic representation 
for bounded harmonic functions h on M as 



h(x)=W[f(A ( )}=E* 



lim/(pr 3 (X t )) 

t >c, 



with / : S* 1 — > E measurable. However, Borbely [6] describes a way to construct 
a family of harmonic functions i/}(r, s) which are rotationally invariant, i.e. inde- 
pendent of a, and therefore cannot be of the above form. As he uses "Perron's 
principle" for the construction, these harmonic functions do not come with an ex- 
plicit representation. Put in the probabilistic framework, we learn from this that 
there must be a way to obtain non-trivial shift-invariant events also in terms of 
the components St and Rt of the Brownian motion. Indeed, as will be seen in 
Theorem 14.21 there exists a non-trivial shift-invariant random variable of the form 

lim \Sf 
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Herein S and R are time- changed versions of S and R and q : K+ — * R is a function 
already constructed by Borbely, whose properties are listed in Lemma 14.11 below. 
This finally leads to additional (to that depending on the component a) harmonic 
functions via the stochastic representation 



h(x) = E x 



g lim I St- 



ilt 



q(r) dr 



where g: 



is a bounded measurable function. 



Lemma 4.1. There is a C°° -function q: 
(1) 



with the following properties: 



q(r) 



sinh(r) 



cosh 2 (r) \ \fh 



— forr<Ti 



(2) For r > Ti the function q satisfies the inequalities 



-3\q\ <q'< 



cosh(r) ' 

(3) There is a T2 > T\ such that 



1 

Vh 



. .Po 40 

< o < . 

" q ~ 2 h 



Mr) 40 

q{r) =—-W) forr - 2 ' 

Proof. See Section l4~2l below and [6]. 

Theorem 4.2. Consider Zt = St — f^* q(u) du as before. The limit 

Zoo ■= lim Z t 

t — >oo 



□ 



Me 



exists. Moreover the law of the random variable (Z^Aoo) = {Z^,Aq) has full 
support IxS 1 and is absolutely continuous with respect to the Lebesgue measure. 

Proof. Let f(r, s) = log g(r, s) and 

9 

m(r, s) 



vW%)HM 

Then a straight-forward calculation shows that 



(4.1) 



dZ, = [ -p(R t ,S t )^(R t )m 2 (R t ,S t )-^q'(R t )m 2 (R t ,S t ^j dt 

(p(Rt,St) - q(Rt)) dt - q(Rt)m(Rt, S t ) dW x t + d W t 2 . 



When R t is sufficiently large, this simplifies as 
dZ t = 



h{R t ) 



40 



(4.2) 



p(Rt, S t )^-{R t )m 2 {R u S t ) - \q'{R t )m 2 (R u S t ) 

4ft I h(R t ) 

q(Rt)m(R t , S t ) dWl + dW 2 . 



dt 



h(R t ) 



For t sufficiently large, we know that \R t — t\ < [3t (for any small (3) and that St 
is positive; furthermore t' r is larger than 2coshrsinhr for r large and s > 0. This 
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together with the fact that the functions p, (h' r /Ah) q', q and 1/Vh are bounded 
allows to conclude that Z is converging. 

To prove that the law of Z M has full support on K, it is sufficient to establish 
convergence of the finite variation part of Z, and that its diffusion coefficient is 
eventually bounded below by a continuous positive deterministic function and above 
by e~P l for some (3 > 0. Indeed, with these properties it is easy to prove that, taking 
any open non empty interval / of K, there is a time t such that after time t, the 
process Z will hit the center of I and then stay in / with positive probability. 

The convergence of the drift has already been established, together with the 
upper bound of the diffusion coefficient. To obtain the lower bound, we use Eq. (|5.8|) 
of Lemma ISTTl below, along with the fact that eventually t/2 < Rt < 2t and St < t, 
to obtain 

for some Ci > 0, where / is defined at the end of Sect. O This immediately yields 
a lower bound for the diffusion coefficient h(Rt)^ 1 ^ 2 m(R t , St) of the form 

((pofr)°/)(2t,t)V 1/a 



ce~ 2t 



(Jio/)(2t,t) 



( Po h) (t/2) 

(for some c > 0) which is a positive continuous function. 

The assertion on the support of the law of (Z m , Aoo) is then a direct consequence 
of the fact that conditioned to (R, Z), the random variable A^ is a non-degenerate 
Gaussian variable. 

We are left to prove that the law of (Zoo, A») is absolutely continuous. Using 
the conditioning argument for Aoo, it is sufficient to prove that the law of Zoo is 
absolutely continuous. To this end we use the estimates in Sect.[SJ It is proven there 
that there exists an increasing sequence of subsets B n of f2 satisfying U n> g B n a = 
such that on B n , for each n > 0, the diffusion process 



U u := (-Rtanu - tan it, Z x 



tan u ) 



has a continuous extension to u S [0, tt/2] and that this extension has bounded 
coefficients with bounded derivatives up to order 2. Consequently, following [Hi] 
Theorem 4.6.5, there exists a C 1 flow tp(ui,v,2,u>) which almost surely realizes a 
diffeomorphism from ]— tanui,oo[ x R to ]— tan«2,oo[ x M if U2 < tt/2 and from 
]— tanui, oo[ x K to its image in R x K if 1*2 = tt/2 (for the first assertion one uses 
the fact that Brownian motion has full support). 

As a consequence, the density p^u of 0^/4 (which exists and is positive on 
]— tan7r/4, oo[ x R) is transported by the flow to 

(0(tt/4, tt/2, c^))- 1 I det J (0(tt/4, tt/2, to)) 1 | 

where J (0(tt/4, tt/2, u>)) 1 is the Jacobian matrix of (0(tt/A,tt/2,uj)) . Integrat- 
ing with respect to uj proves existence of a positive density p^/2 for U v /2- 



P-r/2 



E 



Ptv/4 o (0(tt/4, tt/2, • )) I det J (^(tt/4, tt/2, ■))" 



Finally projecting onto the second coordinate gives the density of Zoo, which proves 
that its law is absolutely continuous. □ 
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4.2. Construction of the function q. As the explicit construction of q is already 
done in [6] we just give a short sketch (following Borbely) how to get a function q 
with the required properties: 

Let a e M+ and T such that p (r)h(r) > 240 and y/h(r) > 80 for r > T . Let 
further T x > T such that p(r, s) = p (r)/2 for r > T x and s > a - 1. For r < Ti 
the function g is defined as 



q(r) 



1 \ ' sinh(r) 



\/hJ cosh 2 (r) 



For r > Ti choose a strictly increasing C°°-extension of — sinh(r)/cosh 2 (r) such 
that q(r) > for r large enough and 



Vh) < q < Vh' 



As p (r)/2 - 40/ft(r) > with 



i' 1 / % 40 , 
Lm (^ 0(r) -%))=° 



r — >oo 



(due to the construction of po, see Sect.QJ there is an r > T\ with q(r) = po(r)/2- 
h(r)/40. Let 



For r>T 2 set 



1 40 
T 2 := inf <{ r > T x : q{r) = -p (r) - 



q( r ) ■= oPo(t-) 



2" uvy /j(r) 

The desired function q is then a smoothened version of this function. 

Borbely shows in 6 , p. 232, that the function q obtained this way meets indeed 
the additionally required properties 

-M<q'< Th and [ Th )<^<- 2 P0- T . 

4.3. Interpretation of the asymptotic behaviour of Brownian motion. We 

conclude this chapter with some explanations how the behaviour of the Brownian 
paths on M should be interpreted geometrically. 

As seen in Sect. [T]and Sect. HI the random variables 

Rt 



lim A* and lim St — q(r) dr 
t-i \ Jo 



serve as non-trivial shift-invariant random variables for X and hence yield non- 
trivial shift-invariant events for X. 

To give a geometric interpretation of the shift-invariant variable 

hm fat - P q{r)dr\ , 
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we have to investigate again the stochastic differential equations for St and R 
dR t = dt + _ dWl 



gh' r + g' r h 











1 










4 9 j 



We have seen in Sect. 2] that the local martingale parts 

,4/i + 4 ff 



M i = J * f^. + £\ dW i and M ?= £ (h(^ + ^)) -/I I J 



of Rt and S* converge almost surely as t — > C- This suggests that the component i?t, 
when observed at times t near 00 (or when starting X near L(oo)), should behave 
like the solution r(t) := ro + t of the deterministic differential equation 

r = 1. 

From the stochastic differential equation above we get 

& = S + — 9 J {Rt ' St l ^^ds + Ml 

Jo g(R t ,St)K(Rt) + h(R t )g' r (Rt,S t ) 

where the local martingale M? converges as t — > 00, and Rt is expected to behave 
like ro + t, when the starting point (?*o, so,ao) of X is close to L(oo). One might 
therefore expect St to behave (for t near to 00) like the solution s(t), starting in sq, 
of the deterministic differential equation 

u ^ . g,(r(f),s) 

1 • J ff (r(i),a)ft{.(r(t))+^(r(t))flJ.(r(t), a )- 

It remains to make rigorous the meaning of "should behave like" . 

Considering the solutions r(t),s(t) of the deterministic differential equations 
above, we note that r so : M + — > K + x K given as T So (t) :— (t,s(t)) with r SD (0) = 
(0, So) is the trajectory of the "drift" vector field 

^ Vd= £h : + gh' r + hg' r !h 

starting in (0, sq) — L(sq). 

As we are going to see below (see Remark the "endpoint" 

r so (+oo) = lim r SQ (t) 

t — >oo 

of all the trajectories is just L(oo) € Soo(M). Furthermore, for every point (r, s) S 
R+ x M, there is exactly one trajectory r so of Vd with r So (r) = (r, s), in other 
words, the union 

S 6R 

defines a foliation of H . Recall that H is one component of H \ L and M — 
(H UL) x g S 1 . Defining a coordinate transformation 

$:I + xI^l + xI 

(4 " >] (r,s) ^ (r,« ) = (# P (r, a), $,(»•,«)), 
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where so is the starting point of the unique trajectory T Sa with r so (r) = (r, s), we 
obtain coordinates for R+ x R where the trajectories r so of Vd are horizontal lines. 

Applying the coordinate transformation $ to the components Rt and St of the 
Brownian motion, 

*(Bt, S t ) = (* r (flt, 5 t ), *.(5t, 5*)), 
we are able to compare the behaviour of the components R t and St with the trajec- 
tories r so of Vd, in other words, with the deterministic solutions r(t) and s(t). The 
component <j> r (i? t , St) obviously equals R t . Yet, knowing that for t — > oo the new 
component <& s (Rt, S t ) possesses a non-trivial limit, would mean that the Brown- 
ian paths (their projection onto (H U L), to be precise) finally approach the point 
L(oo) S Soa(M) along a (limiting) trajectory T so , where s = lim t _ 00 & s (R t ,S t ). 
This would contribute another piece of non-trivial information to the asymptotic 
behaviour of Brownian motion, namely along which trajectory (or more precisely, 
along which surface of rotation r so x S 1 ) a Brownian path finally exits the mani- 
fold M. 

It remains to verify that the above defined component <f> s (Rt, St) indeed has a 
non-trivial limit as t — > oo. As already seen, <3> s (i?t, St) is the starting point of the 
deterministic curve s(t), satisfying the differential equation (|4.3II with s(Rt) = St- 
The solution s(t) is of the form 

s(t)=s + f(r(u),s(u))du 



where / = g' s /(gh' r + hg' r ). In particular, s(t) explicitly depends on s(u) for u <t. 
That is the reason why, when applying Ito's formula to § s (R t ,St), there appear 
first order derivatives of the flow 

$:M + xI^ R + x R, (r,s)^T s (r) 

with respect to the variable s. Estimating these terms does not seem to be trivial, 
nor to provide good estimates in order to establish convergence of & s (Rt,St) as 
t — > oo. 

A possibility to circumvent this problem is to find a vector field V on T(R + x R) 
of the form d/dr + f(r)d/ds whose trajectories also foliate H and are not "far 
off" the trajectories r So of Vd - in particular the trajectories of V have to exit M 
through the point L{oo) € Soo(M) as well. 

As seen in Sect. SI we have 

I 



9's 

P 



gh' r + hg' r 



< 



P ■ 



l + h 



Furthermore, for r > T2 the function q(r) is defined as po/2 — 40//i, in particular 
q(r) does not differ much from the function p(r,s) which equals po/2 for r and s 
large. Hence q(r) is a good approximation for g s /(gh' r + hg' r ) for r large, and is 
independent of the variable s. 

We therefore consider the vector field 

Starting in (0, s ) € R+ x R the trajectories C so of V have the form 

C SQ (t)= [t, so+ / q{u)du 
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Figure 4. Effect of the coordinate transformation <f> 



As we are going to see below, we also have limt_ s . 0O C' So (t) = L(po), see Remark l4.3[ 
and the union 

U Coo 

s 6R 

forms a foliation of H. 

For (r, s) € R+ x R there is exactly one trajectory C Sa of with C So (r) = s. Its 
starting point so can be computed as sq — s — f Q q(u) du. We can therefore define 
a coordinate transformation 

$ : R + x M -> R + x R, 



(4.7) 



(r,s) 



r, s 



As seen in Figure [H the trajectories C So of V are horizontal lines in the new 
coordinate system. 

In the changed coordinate system the components Rt and St of X then look like 
(4.8) ^(R u S t )=(R u St- J \{u)du) = (® r (R t ,S t ),<S> s {R t ,St)). 



As we have proven in Sect. 0] 
liin St — 



q(u) du 



lim$ s (i? t ,S t ) 



exists and is a non-trivial shift-invariant random variable. Therefore the non- 
triviality of lim^oo <I> s (i? t , St) allows to distinguish Brownian paths when examin- 
ing along which of the trajectories C So of V, or more precisely along which surface of 
rotation C S() x S 1 , the path eventually exits the manifold. This gives the geometric 
significance of the limit variable 



lim 

t — >oo 



Si 



Rt 



q(u) du 



It finally remains to complete the section with the proof that the trajectories of 
the vector field V, as well as the trajectories of the vector field Vd, exit the manifold 
M in the point L(oo). This is done in the final remark. 
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Remark 4.3. We have: 



lim C so (t) = L(oo) and lim r so (t) = L(oo) for every s € K. 

t — >oo t — >oo 

Proof. It is enough to show that the "s-component" of each trajectory C So , T So 
resp., converges to oo with t — > oo. The s-componcnt of C sa is 



s(t) =s + / q(r)dr, 



the s-component of r so 

s(t) = s Q + 
Since for t > T2 we have 



ffa(r,s(r)) 



st^sWl^rH +/i(r)fl^.(r,s(r)) 



dr. 



( \ 1 ( \ 40 



it follows immediately that lirnt_,oo f q{r) dr — 00 because J °° ^4pr dr < 00 and 



lim 

i — *oo 



Po(r) dr = 00, 



due to Lemma [2TT1 property ([3]). 

For the second term we notice that so + L* /(f, s(r)) ^ r with / = g's/(gh' r + hg' r ) 
is nondecreasing as the integrand is positive. Moreover we have seen above and in 
the foregoing sections that 



As 



gh' r + hg' r 
p(r, s) 









p 


< 


p -l+h 



1 



1 + h(r) 
t 



dr < 



l + h(r) 



dr < 00 



it suffices to show that lim^oo J Q p{r 1 s(r)) dr = 00. This is true as s(r) > sq for all 
r < t and therefore for r large enough we have p(r, s(r)) = po(r)/2. The claimed 
result then follows exactly as above. □ 



5. The Poisson boundary of M 

In this section we prove that any shift-invariant event for Brownian motion in M 
is measurable with respect to the random variable (Z^,A^) constructed in Sect.|U 
As a consequence this allows a complete characterization of the Poisson boundary 
of M. 

We perform the change of variable \& : (r, s, a) 1— ► (r, z, a) with z — s — J q(u) du. 
Further let t(r, s) = log g(r, s) and 

( \ 2 
m(r,s) = — ^^^^=^^^^= . 
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In the new coordinates, the three components Xt(x) = (Rt(x), Zt(x), At(x)) = 
(Rt, Zt,At) of time-changed Brownian motion satisfy 

dR t = dt + m(R t ,S t )dW t \ 

dZ t = ('-p(^,5t)^(H t )m 2 ( J R t ,5 t ) + -^ / (^)m 2 (5 t ,^)^ dt 

(5.1) + ( p (R u S t ) - q(Rt)) dt - q(R t )m(R t ,S t ) dW x t + "H^M dW ^ 

\Jh(Rt) 

d A t = ^£Ldwf. 

\j g(Rt,S t ) 

In case when R t exceeds the positive constant Ti defined in Lemma |4. 11 the equa- 
tions simplify to 

dR t = dt + m{R t ,St)dWt, 

dZ t = m 2 (R u S t ) (-^o(i? t )^(i? t ) - dt + ^- ) dt 



( 5 - 2 ) + m(R t ,St) \ -q(Rt) dWl + = dW? | , 

'h(Rt) 



dA t = ^feg dW?. 
y g(Rt,S t ) 

In the subsequent estimates we always assume r > T 2 and s > 0. Since p is 
nonincrcasing and p h is increasing and strictly convex, we have (pgh)' > which 
yields 

h! 

0>Po > -Poj 

and consequently 

(5.3) \p' (r)\ < 2po(r), |<z'(r)|<PoW- 
It is easy to see that 

(5.4) \(h- 1/2 )'\ < h- 1 ' 2 . 

We want to estimate m, m' r and m' s . We know that t' s = pht' r . Recall that for 
So > 0, ro > 0, the curve s <— > (7(^0, So, s), s) is the integral curve to the vector field 
d s — ph(r,s)d r satisfying 7(ro, so, s o) = J"o- Recall further that /(ro,so) is defined 
by 

f(r , s ) = 7 (r , s , -l{f{r , s ))) , 
see Eq. (|2.10l) . and that by definition of the metric (7, we have 

ff(r 0) so) = 50 (f(ro, *o)) where 5o (r) = d 2 e dl sinh2 r 
for some di, cfe > 0, see Eq. (|2.11[) . Letting 

t (r) = log 50 M, 
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the function t(r, s) = log g(r, s) = to (/(»", s)) satisfies 

(5.5) t(r ,s ) = di sinh 2 (/(r ,s )) + logd 2 , 
which yields 

(5.6) \t' r \ <2d 1 {ho f)\f r \ and \t' s \ = ph\t' r \ < 2d lP h(h o f)\f' r \. 

Lemma 5.1. There exists two constants C\,Ci > 1 such that on the set {(r,s) £ 
[3,oo[x [0,oo[}, 

(5.7) Mll<f r < Cl ^ll t 

p h p h 

and 

(5.8) 6 2 ho} — <t r <C 2 hof — . 

p a h p a h 



Proof. It suffices to verify (|5.7p . Indeed, assuming that (|5.7p is true, then from the 
equality t' r = (t' o /)# and the bound h(r)/C" < t' (r) < C'h(r) for some C > 1 
and all r > 3, we obtain (|5.8|) . 

To establish (15. 7|) . let (r ,s ) € [3, oo[ x [0,oo[ and denote by 7- the derivative 
of 7 with respect to the i-th variable, i = 1, 2, 3. We then have 

(5.9) 73(7-0, s , s) = -(p/i) (7(r , s 0) s), s) 
and 

(5.10) 7(r , s , -£(f(r , s ))) = /(r , s ). 
Differentiating (|5.10[) yields 

/r( r O,So) 

= 7i(n), s , -£(f(r , s ))) - 73( r o, s o, -t(f{ra, s )))£'(f(r , s ))fr(r , s ) 

= 7l (r , so, -*(/(r„, s )) + (ph) (-£(f(r , s )), /(r , s )) ^(/(r , «o))/;(n>, s ) 

= 7 1 (r ,s ,-^(/(r ,so)), 

since ph(r, s) = x(r, s){p h)(r) and x (/(r"o, s ), -£(f(r , s ))) = (see [5] p. 229). 
Consequently it is sufficient to prove that 

(Poh) o f , , , , ^ (Po^) o/ , s 

(5.11) - — (ro,so) < 7i(?"o,so,— K(/(r ,s )) < 6 - — (r ,s ). 

p h p h 

To this end, we differentiate (|5.9[) with respect to r and obtain 

(5.12) 7i3(n>,so,s) = -(ph)' r (■y(r ,s ,s),s)'y[(ro,so,s). 
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Solving the last equation with initial condition j[(ro, so, so) = 1 (which comes from 
7(7-0, s , s ) = r ) gives 

7i (rQ,s ,-£(f(r ,s ))) 
= exp[/ {ph)' r (j(r ,so,s),s)ds ] 

\J-e(f{r ,s„)) J 

= exp( f x(l(rQ,so,s),s){p hy('y(r ,so,s))ds) 

\J-e(f(r , S0 )) J 

x exp ( f Xr(l( r o,so,s),s)(p h)(~/(r ,s a ,s))ds) 
\J-e(f(r .s )) J 

( f S0 h< f , , bofey( 7 (ro,so,s)) \ 

= exp / ph{-y{r ,s ,s),s)- — — — : w ds ] 

\J-£(f (r ,s )) [Poh)h{ro,so,s)) J 

x exp ( f Xr(l( r o,s ,s),s)(p h)(j(r ,s ,s))ds\ 

( [ So f , (poh)'(>y(ro, So , s )) \ 
= exp / -73 r ,s ,s) 77 ds 

x exp ( [ Xr(j( s o,r 0l s),s){poh){~/{ro,S( h s))ds) 

where for the last equality we used (|5.9|) . We thus get 
7i (nb«o,-^(/0"o,so))) 

= exp ( f -— (log(p h)(j(r , s , s))) ds ] 

V-<(/(r 0) »o)) rfS / 



x exp / x'ri'y{ r o,so,s),s)(poh)(-f(r ,s ,s))ds ) 

IH,L () ^exp( / Xr(7( 7 'o, so, s),s)p Q h(~f(r , s , s))ds 



MM \./-*(/(m,.o)) 

Poh(f(r ,s )) 



exp / ^ / (s + ^( 7 (ro,so,s)))(£ / po/i)(7(ro,so,s))ds 



MM V-^-^/^o.so)) 
Poh(f(r Q ,s )) 



exp / ^'(s + £( 7 (ro,so,s)))eds 



MM \/-<(/(n„«o)) 

since for s e [-£{f(ro, s )), s ], 7 (r ,s ,s) > 7 (r ,so,s ) = r > 3 and f(r) 
s/(poh(r)) for r > 3 (see [6] p. 229 and Sect. [2]). Hence we have 

7i Ms ,-£(/Ms ))) 

M(/Ms )) ( r-W™))** \ 
■exp / £ (s + ^( 7 (r ,So,s)))eds 



MM U_« (/(ro , So) ) 
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where we used the fact that £'(s') = when s' > 4, and 

^(s + £( 7 (r , s , a))) = 1 - e'ph(j(r , s , s)) > 3/4 

(0 (2.15) and e < 1/4). 

Since £' is nonnegative and bounded by 1/2, we have 

/ f-e(f(r ,a ))+8 \ 

1 < cxp / £(s + £("f(r ,s ,s)))eds < exp(4e) 

\J-i(f(ro,s )) J 

which finally gives 

(5.13) M^(r , So ) <ii(r ,s ,-£(f(r , So ))) < exp(4 £ )^^(r , So ). 

poh poh 

This is the desired result. □ 
Lemma 5.2. There exist constants C\, C% > such that 

(5.14) \a<c u poh 

and 

(5.15) | C |< Cbfco /(^Z) a . 

Proof. Assume that (|5.14l) is true. From the identity 
(5-16) C = (*o°/)(/;) 2 + (*o°/)^ 

along with (|5J)) and the fact that < t' Q < 2d 1 h, < t' Q ' < 4dih, we obtain ([535)) . 
To establish (|5.14|) we first note that 

/r( r o,s ) = 7i( r o,«o,-^(/(ro,so)) 

and 

7iOo,so,-^(/(r ,s ))) 
(5.17) = — — exp / f(s + €(7(r 0) s ,s)))eds . 



Poh(r ) \J-e(f{ra,s )) 

Hence, we have 
(5.18) 

/"r( r O,So) 

_ ft-./.)' (/ft,, «,)) KM ( r «-)>« ao _ A 



MM V-WCm,*))) 

/ — tto? \ exp / f (s+£(7(ro,s ,s)))£ds 

+ /;(r , 8 ){f(f(r ,8o))K(r , 8 )?(0)) 

»_^(/(r , JSo ))+8 

+ /r( r o-s ) / ^"(s + ^(7(r ,so,s)))^(7(ro,so,s))7 1 (r ,so,s)eds. 
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Using < (poh)' < 2poh, along with the boundedness of the integral in the expo- 
nential, the bound for f' r , the boundedness of £" and £' ([6] (2.15)), and the fact 
that 

< j'i(r Q ,s ,s) 
= 6XP (/ <yph ^ r ^^ ro,So,u ^ u ^ du ^) 

< exp / (ph)' r (~f(ro,so,u),u)du\ (since (ph)' r > 0, see p. 229) 

\J-e(f(r , So )) J 

= j'i(ro,so,-£(f(r ,s Q )) = f' r {r ,s ), 
we get the wanted bound. □ 
Lemma 5.3. There exist f3 > 1 and vq > such that for all r > ro and s > 0, 
(5.19) f(r,s)>(p Q hf(r). 

Proof. Since f(r, s) > f{r, 0) for s > 0, it is sufficient to establish (|5.19p for s = 0. 
Let r be such that i(ro) > 4. For r > r , define /i(r) as 

7(r,0,-*(r) + 4) = /i(r). 

Since —£(f(r,0)) < —£(r) + 4, we have f\{r) < f(r, 0). Hence it is sufficient to 
establish (|5 . 19[) with fi(r) in place of f(0, r). Writing 

fh(r) A 

-£(r)+4= / — ~/- 1 (r,Q,-){u)du 

Jr du 

we obtain 

f (r) - 4 = / — — r- — — du. 

ph(u,1 {r,0,-){u)) 

But when u G [r, /i (?")], we have 

7~V, 0, ■)(«) > 7 _1 ( r 7 0, 0(A(r))) = -*(r) + 4, 

hence 

£(u) + 7- x (r, 0, > £(it) - £{r) + 4 > 4, 

and this implies 

7~ V, 0, ■)(«)) = + 7~ 1 (r, 0, •)(«)) = ~ 

Thus, since p/i(r, s) = x{ r i s )Po( r )M r )j we obtain 

2 

p h(u) 

yields 

/(r) 



/■AW o 

(5.20) £( r )-4= / - du. 

J r 

Differentiating with respect to r yields 

2/l(r) 2 



PoHfi(r)) Poh(r)' 
which combined with £' = e/(poh) gives 

(521) (2 + s)p h(f 1 (r)) 

(J - 21) 2p h(r) ■ 

From the convexity of p^h it follows 

Poh(fx(r)) > (p h)'{r) (h(r) - r). 
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By Eq. (15.20|) . taking into account that poh is nondecreasing, we derive the inequal- 
ity fi(r) - r > (e(r)-4)p h{r). Choosing A € ]0, 1[ such that f3 := A (2 + e)/2 > 1, 
this implies fi(r) — r > A/i(r) for r sufficiently large. We get 

P*h{h{r))>(p Q h)'{r)\h{r) 

and with ([5T2"Tj) 

(5-22) — — > /3 — — . 

fi(r) poh(r) 

Integrating (|5.22[) from tq sufficiently large and noting that /i(j~o) > po^( r o)) we 
finally arrive at the desired result. □ 

We are now in a position to estimate m' r and m' s . Since and its derivatives 
are bounded, we can estimate \m' r \ by C \d r (t' r )~ 1 / 2 \. Recall that t(r, s) = t (f(r, s)) 
which gives 

(0- 1/2 = (*o°/)- 1/2 (/;)- 1/2 

and 

dr(t' r r 1/2 = o /r 3/2 (^ o /)(/;) 1/2 - |c*& ° zr 1/2 (/;r 3/2 (/;;). 

The last equation, along with ([577j) and f5~T4")) . and (1/C)h < t , t' Q , t'^ < Ch for 
some C > 1 and all r sufficiently large (recall that to(r) = d\ sinh 2 r + logc^), gives 

(5.23) \m' r \ <C(p °/) 1/2 (Po^r 1/2 . 

Alternatively, using Eq. (|5 . 23[) together with Lemma [5.11 and estimating mT 1 by 
(4) 1 / 2 , we obtain 

(5.24) !<! < C {P ° k) ° f . 

m (p h) 

Similarly to m' r , we estimate \m' s \ with C \d s (t' r )~^ 2 \. Since = (ph)t' r , we have 

= -i(0" a/2 a rfipo K 



2' 

-3/2 4 



Recall that |(po^)'| < 2poft-- The first term on the right can be bounded in absolute 
value by 

cxwr 1/2 «M) ° /r 1/2 (M) 3/2 

which is smaller than C(po ° fY^iPoh) 1 / 2 since 

poft < p h < 



By means of Lemma 15.11 and Lemma 15.21 the second term on the right can be 
bounded in absolute value by 

COW) 1/2 (M) 1/2 • 
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Consequently we obtain 

(5.25) \d s (t' r y 1/2 \ < C(p ° f) 1/2 (p h) 1/2 . 

The following lemma is a consequence of Lemma \5 .'61 and Eqs. (|5.23p . (|5.25[) . 
Lemma 5.4. There exists n > such that for every r > r\ and s > 0, 

(5.26) K| < C( Po o /)Va ( Po fc)-Va j | TO > | < C ( Po o (p^h) 1 ' 2 , 

(5.27) ^ C ^ — rr < C([p h) o f). 

m (Poh) m 

We recall that 

(5.28) C- 1 {hof)-^ u y' 7 ' <m<C(hof)-^- iPah)1/2 



((p Q h) o /)i/a - - v " ((poM Q /)1/2 

for some C > 1. Let us briefly explain how bounds for higher order derivatives can 
be established. Differentiating expression (|5 . 18[) with respect to r yields, for s > 
and r sufficiently large, 

'(Poh) o p 



(5.29) <C nj! 

(We remark that differentiation of each term in (|5.18p amounts to multiplication 
by an expression smaller than C((poh) o f)/(poh) in absolute value). This easily 
yields, for s > and r sufficiently large, 

(5-30) \4%\<C(hof)(^^iy 
and 

(5.31) \rn';A<C(hof)-^(^^y\ 

Now exploiting the fact that t' s = pht' r , a straightforward calculation shows (for 
s > and r sufficiently large): 

(5.32) \ t W\<c(poh)(hof)^^r a 

(5.33) \tit\<C( P0 hnhof)^^ 

(5.34) K r | < C( Po h)(h o /)-Va (l^/^ 2 , 

(5.35) KJ < C(po/i) a (fto/)-Va (1^T) 3/2 , 

(5.36) K r m| <C{p h)-\nof), 

(5.37) K>| < C(p o /), 

(5.38) \m'; s m\<C(p h)(p °f). 

Proposition 5.5. For any a € ]0, 2(0+1) [ there exists a constant r a > suc/i iftai 
/or all r > r a and s > 0, 

(5.39) |m|, K|, |m' s |, |mm£.|, |m<,|, |mm£| < h~ a (r). 
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Proof. All the estimates are either immediate or easy consequences of the ones 
for \m' s \ and |mm" s |. On the other hand, since the proof is the same for these 
two functions, we only establish the estimate for |mm" s |. Note that the final a 
should lie in ]0, 2(8+1) [ wnere /3 > 1 is given by Lemma [5.31 due to the estimate 

K| < (po/0 1/2 0W)V2. 

Fix a £ ]0, 4qpj[. It is sufficient to prove that h a (r)(poh)(r)(p ° f)(r,s) is 
bounded for r > ro and s > 0. Then the claimed result is obtained by picking a 
smaller a. 

From Lemma 15.31 and po < 1/r, we have 

h a (r)(p Q h)(r)( Po of)(r, S ) < h a {r){p h) 1 -^ \r). 

If /i a (r)(pQ/i) 1_/3 (r) < 1, we are done. Otherwise we have (p h)(r) < /i"^ 1 ) (r) 
which yields 

p (r) < W^- 1 )- 1 ^). 
In this case, since f(r, s) > r and po is nonincreasing, 

h a (r)( Po h)(r)(p a of)(r, S ) < h a (r)h a 'V > - 1 )(r)h a 'W- 1 '>- 1 (r). 
We finally get 

h a (r)(p h)(r){p o /)(r, s) < 1 V (ftW/5+l)-(/3-l)]/(/3-l)( r ) 



since a £ ]0, the right hand side is clearly bounded. □ 

5.1. Equation for the time- reversed process. Consider 

X t =X t (x) = {R t ,Z u A t ) 

the time-changed Brownian motion started at x, satisfying Eq. (|5.1|) . 

Let t — tanu and define X u = (R u , Z u , A u ) = X ianui < u < n/2. Furthermore 
denote 

m(r, z) — mo (r, z, a) = m(r, s), 
p(r,z) =por 1 (r,z,o) = p(r,s), 
g(r, z)= go * _1 (r, z, a) = g(r, s). 

From (|5.1|) we get 



dR u = — \—du H —m(R u ,Z u ) dW^, 

cos z u cos 14 



(5.40) 



dZ„ = — \— ( -p{R u ,Z u )^(R u )m 2 (R Ul Z u ) - \q'(R u )m 2 (R u ,Z u ) ] 
cos z u \ 4n, 2 / 

1 _ u 

{p{R Ul Z u ) - q(R u ))du 



cos 2 u 



1 ■ A)^, Z u ) dti£ + J- dWl 



COS U COS 14 1 , r r, \ 

Jh(R u ) 

g{R u ,Z u ) 
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for some three-dimensional Brownian motion (W 1 , W 2 , W 3 ). Letting Y u = R u 
tanu and X u = (Y U7 Z U , A u ) — (X*, X*, X 3 ), we can write 

(5.41) dXl =ai(u,X u )dWi t +V(u,X u )du, j = 1,2,3, 

where b 1 = b 3 = 0, 

1 h' 

b 2 (u, x) = 7i — pix 1 + tanw, x 2 )—^-{x 1 + tanu) m 2 (x 1 + tanu, x 2 ) 

v ' cos 2 u v ' 4/i v ' v ' 

■ q ' (x + tan it) fh 2 {x 1 + tanu, x 2 ) 



9 o 

cos u 2 

H ^ — {pix 1 + tan it, x 2 ) — qix 1 + tan it)) 

cos 2 u 

40 

(the last line equals - — — — — r- for x 1 + tan it large), 

cos 2 u h{x L + tan it) 

a}(u, x) = mix 1 + tan it, x 2 ), a\ = a\ = 0, 

iv ' cosu 6 

cxAu, x) — qix 1 + tan u) mix 1 + tan it, x 2 ), 

cos u 

1 mix 1 + tan it, a; 2 ) 



a 2 {u,x) 



J 2 ~3 



<J 3 {u,x) 



cosu yjhix 1 + tanu) 
1 rf^x 1 + tantt, a; 2 ) 



cosu ^g(x x + tanu,x 2 ) 

Adopting Proposition 15.51 and defining a{ir/2,x) = and b{%/2,x) = 0, we obtain 
easily a C 2 extension of a and 6 on D = {(u, x) € [0, n/2] x R 3 , x 1 > — tan it} with 
vanishing derivatives at u = tt/2. 

Let K be a compact subset of D. Further let a K , b K be two C 2 maps defined 
on [0,7r/2] x IR 3 which are bounded together with their derivatives bounded up to 
order 2, and which coincide on K with a and b. We fix a £ ]0, 7r/2[ and let X K be 
defined by 

(5.42) dX* = o K {u,X u x )dW u + b K {u,X u x )du with If = X a . 

For the diffusion (X„ )„ e [ Q . ]7r /2] conditions (HI), (H2) of [55] are satisfied, and we 
may conclude with corollaire 2.4 of that the time-reversed process 

{Xu = X^ 2 _u)u6[0,jr/2-oi] 

satisfies 

(5.43) dX* = <7 K (u, X* ) cW u + Xf ) du with Xf = Xf /2 , 

where W is a Brownian motion independent of Xg and where the coefficients are 
given by 

a K {u,x) = cr K (7r/2-u,x), 

(6 K (it, x)Y = -(6 K (7r/2 - u, a;))' + (p A '(tt/2 - u, a;))" 1 ^-({ a K )^ P K ){-K 12 - u, x); 

ox 1 

here we have {a K ) l i = ELiCOU "^)! and K x ) is the density of X* . 
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5.2. Invariant events. 

Proposition 5.6. For every x = (r, z, a) 6 M with r > 0, up to negligible sets, the 
pre-image of J^inv under the path map $ to the Brownian motion X = (R,Z,A), 
starting from x, is contained in the a-algebra a{X^/2) generated by X^/2- 

Proof. Fix a € ]0, 7r/2[. Since ^i nv C <^oo, the cr-algebra $ _1 (^i n v) is clearly 
included in P <r(X s , s > <) = $ _1 (^ 00 ). On the other hand, we have 
t>o 

P| a(X s , s > t) = p| a(X v ,v>u) 

t>0 ue]a,7r/2[ 

= P c(X„, w > u) up to negligible sets 

u£]a,7r/2[ 
= P cr(X„, t) < u). 

ue]0,7r/2-a[ 

For n > 1, consider the compact subset -K„ of D, 

K n = \(u,x) G [0,7i72l x M 3 , a; 1 > -tanu + -, \x l \ <n, i = 1,2,3 
and let 

B n = {cj, e K n , Vu € [a,7r/2]}. 

Since X has a.s. continuous paths with positive first component, we have for any 
n > 1 

(5.44) |J B m ^ n . 

m > n 

Now define X n = X K ™ and X" = X K ™ . For 

A e P <t{X v , v < u), 

«e]0,7r/2— a[ 

we find from Eq. (|5.44[) 

(J AnB„. 

n>l 

Denoting by $ (resp. <£„) the path map corresponding to X (resp. X n ) and by 
the cr-field of terminal events for continuous paths [a, 7r/2] — > D, there exists 
C G &tt/2 such that 

A = $- 1 (C'). 
Then on B n , we have X a =' X n , and hence 

(5.45) AnB n $; 1 (C)ns n 

where 

^(CJe P 

ue]a,7r/2] 
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On the other hand, since X n is the strong solution of a stochastic differential equa- 
tion driven by a Brownian motion W n independent of Xq (see [23] corollaire 2.4), 
we have 

p| *(X2, v>u)= f] a(X:, v < u) 

«e]a,7r/2[ ue]0,7r/2-a[ 

= P| (a(XZ)Va(W:, «<«)) 

ue]0,7r/2[ 

= <j(X'n ) V Pi cr(WC, u < it) (up to negligible sets) 

«S]0,7r/2[ 

= a(XZ) = a(X: /2 ). 
As a consequence, there exists a Borel subset E n of R 3 such that 

®-\C) = {X^ /2 G E n }. 
This yields, along with (|5.45|) . 

(5.46) A n B n ^ {x™ /2 g £„} n B n ^ {X n/2 e £„} n S„ c {i x/2 e E n }. 

The inclusion AnB„C {X w / 2 S S ra } together with (|5.44|) yields for all n > 1, 

A C < G |J £ m I a.s. 

This is true for all n > 1, so letting _E = limsup_E n , we get 

n — >oo 

A C {X„/ 2 G -E} a.s. 

Now let us establish the other inclusion. If w G {X^/ 2 G E 1 } then choose n such that 

(u,X u (uj)) G if n f° r au u G [of>7r/2]- Then choose to > n so that X„/2(u>) G -E m - 
We clearly have w G B m , so that wein B m by (|5.46[) . Finally we proved that up 
to a negligible set, 

A = {X n/2 eE}. 

This shows that 

f]a(X s , s > t) C a(X* /2 ) 
t>o 

up to negligible sets. □ 
Let Y t = Rt — t and = lim Y t . We proved that the invariant paths for the 

t — ¥00 

process X t = (Rt, Z t , A t ) are measurable with respect to (Yoo, Z^, The last 

step consists in eliminating the first variable. 

Theorem 5.7. Let xq = (ro,zo,ao) G M. TTie invariant paths for the process 
Xt(xo) = (Rt(xo), Z t (xo), A t (xo)) are measurable with respect to (Zoo^xq), A 00 (xq))- 

Proof. Let xq = (ro,zo,ao) G M and S Xo the sphere of radius 1 centered at Xo- 
Define a time-changed Brownian motion X(xq) (time-changed in the sense of a 
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solution to (|5.1|) ) on a product probability space as follows: X'(xo)(u>i) is a time- 
changed Brownian motion started at xq, and for any x G S Xo , X" (x){u>2) is a time- 
changed Brownian motion started at x. Letting r = mi{t > 0, X' t (xo) £ S^}, we 
define 



*t(">i)_ if*<r(wi), 
X t "_ T (X>i))(u, 2 ) ift>r( Wl ). 



Given _B e ^i nv , we need to prove that there exists a Borel subset E' XQ of M x 5 1 
such that P = Pi (g) P 2 -a.s., 

(5.47) {X(x ) eB} = {(ZUxq), 4»(a:o)) G E' Xq }. 

For every x G M, let E x be a Borel subset of K x K x S 1 such that 

{x(x) e B} a ^ {(y^ (x), Zoo (x), A^x)) e e x }, 

where X{x) is a time-changed Brownian motion on M started at x, with coordinates 
(yt(x) + t, Z t (x) , A t {x))) . Since B G ^nv, we have 

{X(x Q )eB} = {X T+ .(x )eB}, 

and hence 

X(x ) S B if and only if X"(X' T ) G 5. 
But by Lemma P-a.s., 

{X(x ) G B} = {(Foo(a;o),^oo(^o),Ioo(a;o)) G £^ }, 
and for all a; G S x „, P2-a.s., 

{X"(x) eB}= {(Y^(x),Z'^(x),A^(x)) G E x }, 
On the other hand, it is easy to see that 

KooteoXwx.ato) = Y^{X' T {ux)){^) - r(wi). 
Let Py, be the law of X^(ieo). A consequence of the equalities above is that for 
P^, -almost all x, we have P( • \X T (xo) = x) a.s., 

{(YM,Z'M,A^(x)) G E x , X T (x )=x} 
= {X"(x)(tu 2 ) G B, X T (x )=x} 

(5.48) = {X(z ) G B, X T (x Q ) = x} 

= {(Yoo(xq), Z oo (x ), Aoo(xo)) G Bx 0) X t (xq) — x} 

= - r(wi), ^(i)^),^,^)^)) G £? X0 , X r (x ) = x}. 

On the other hand, we know that the law of r given X T (xo) and 0J2 has a positive 
density on ]0, oo[ and depends only on X t (xq). Hence equality of the second and 
the last term in (|5.48[1 implies that for Pj, , s-alniost all x, P( • \X' t {xq) = x) a.s., 
if 

(Y^x^-rioj^Z'^x^A^x)^)) G 

then i?a; contains almost all (y', Z'J (x)(ui2), A'^- l {x){ui2)) such that y' < Y£ (x)(u)2)- 
Let 

E' = {(2, a)eMxS' 1 |Elj/eR such that (y, z, a) G B^for almost all y' < y}. 
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We proved that for Py, , r almost all x, Pi( • \X' (xo) = x) ® P2 a.s., 

{X"(x){co 2 ) G 5, X T (x ) = x} 

C {(^(x)( W2 ),2^(x)( W2 )) G X T (x ) = x] 
= {(Zoofeo)^), A oo (x )(w)) G £^ , X T (x ) = x}. 

Since 

{X{x ){u x ,u 2 ) G S, X T (x )M = x} = {X"(x)(w 2 ) G fl, X T (x )(wi) = x}, 
we get almost surely, 

(5.49) {X(x ) £B}C {(Z oo (x ),ioo(x )) G ££„}. 
To establish the other inclusion, let us define, for (z, a) G -E^ ' 

„ ^ _ |sup{y G K, for almost all y' < y, (y' , z, a) G £^ Q }, 
1 —00 if the set is empty. 

We know that on {X(xo) G B}, for Pj, ^.-almost all x, P( • |X^.(xq) = x) a.s., 

/(z oo (x ),i oo (x )) eiuH 

Recall that Y 00 (xq)(uji,W2) = Y^(X^(x )(a;i))(a;2) - r(wi), and that 

{(■£00(^0X^1,^2), Ax>(zo)(wi,w 2 )) G E' XQ , X T (x )(u)i, w 2 ) = x} 
= {(^(x)(c, 2 ),i'^(x)M) G X;(x )(wi) = 1}. 
By the positiveness of the density of the conditional law of r, for Pg, , .-almost 
all x and w 2 G fi 2 such that 

(^(x)(c 2 ),i'^(x)(^))e^ o , 

we find 

Pi {r > - /(^(x)( W2 ),A£,(x)M) I X T (x ) = x} > 0. 

This yields 

P 1 {^( a: )(a; 2 )-r,tWh),^h))e£ IO | X T (x ) = x} > 

which implies by (|5.48p that 

Pi [X"(x)(w 2 ) G B| X T (x ) - x} > 0. 

Taking into account that above probability takes its values in {0, 1}, it must be 
equal to 1. Consequently for Pg, -almost all x, P( • \X t (xq) — x)-almost surely, 

{^(x)M,i'^(x)M) G E' Xo , X' T {x Q ) = x} C {X"(x) G B, X' T (x Q ) = x} 
which rewrites as 

{Z oo (x ),l oo (x )) G E' XQ , X T (x ) = x} C {X(x ) G £, X r (x ) = x}. 
We get, P-almost surely, 

(5.50) {^oc(x ),i 00 (xo)) G ££j C {X(xo) G £}. 
With dEjO]) and (OP} we finally obtain (|CTf)l : 

{Z oo (x ),2 oo (x )) G ££J = {X(x ) G B}, P-a.s. 
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This achieves the proof. □ 

Putting together Theorems 14. 21 and 15 . 71 we are now able to state our main result, 
giving a complete characterization of the Poisson boundary of M. 

Theorem 5.8. Let ^(lx S 1 ) be the set of bounded measurable functions on IRx 5 1 , 
endowed with the equivalence relation f\ ~ fi if f\ = fi almost everywhere with 
respect to the Lebesgue measure. 
The map 

(&(R x S 1 )/ ~) — ► J%(M) 

f^(x»E[f(Z ( (x),A ( (x))}) 

is one to one. More precisely, the inverse map is given as follows. For x G M, let- 
ting K(x, •, •) be the density of (Z^(x), A^(x)) with respect to the Lebesgue measure 
mix S , for all h G J%{M) there exists a unique f G ^(R x S 1 )/ ~ such that 

Vx£ M, h(x) = / K(x, z,a)f(z,a) dzda. 

JRxSi 

Moreover, for all x G M, the kernel K{x, •, ■) is almost everywhere strictly positive. 
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